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SECTION  I 
INTRODUCTION 


Most  solid  propellant  rocket  motors  exhibit  some  degree  of  combustion 
instability,  which  is  characterized  by  chamber  pressure  oscillations.  The 
hot  combustion  gasses  in  the  combustion  cavity  can  oscillate  in  various 
natural  acoustic  modes  much  the  same  way  that  the  column  of  air  in  an  organ 
pipe  resonates.  Special  pressure  transducers,  that  are  designed  to  measure 
the  alternating  component  of  the  chamber  pressure,  are  used  to  measure  the 
unstable  pressure  oscillations.  The  oscillations  are  considered  to  be  un- 
stable because  a small  perturbation  can  excite  a particular  mode  which  in 
turn  increases  in  amplitude  in  an  unstable  (eat  envelope)  fashion  until  some 
limiting  amplitude  is  reached.  As  burning  in  the  motor  continues,  the  condi- 
tions required  to  sustain  oscillation  in  a particular  mode  change  and  the 
mode  typically  dies  away  before  the  end  of  motor  operation.  Motors  that 
exhibit  unstable  pressure  oscillations  in  more  than  one  acoustic  mode  during 
motor  operation  time  are  common. 

In  the  past,  unstable  acoustic  pressure  oscillations  in  upper  stage 
motors  on  certain  ballistic  missiles  have  produced  relatively  high  amplitude 
vibration  levels  on  the  motor  case  and  attached  components.  Vibration  levels 
as  high  as  300  g's  have  been  measured  during  upper  stage  motor  operation. 

This  past  experience  with  solid  rocket  motors  has  been  cause  for  the  concern 
with  possible  acoustic  instabilities  in  the  Space  Shuttle  Solid  Rocket  Booster 
(SRB)  motors.  The  objective  of  the  work  covered  in  this  final  report  is  to 
analyze  the  Space  Shuttle  vehicle  to  determine  structural  response  to  possible 
acoustic  combustion  instability  in  the  solid  rocket  boosters.  The  analysis 
is  to  provide  an  estimate  of  the  forces  to  be  expected  at  the  attachment  points 
between  the  Solid  Rocket  Boosters  and  the  External  Tank  (ET) . 

Work  to  define  the  likelihood  of  any  particular  acoustic  mode  being 
unstable,  to  define  the  natural  mode  shapes,  and  to  estimate  limiting  ampli- 
tudes of  unstable  modes  is  being  carried  on  at  the  Naval  Weapons  Center  (NWC) 
at  China  Lake,  California  under  the  direction  of  Dr.  Ron  Derr.  For  the  work 
reported  herein,  a mode  was  assumed  to  be  unstable  and  the  response  was  calcu- 
lated for  a normalized  (1.0  psi  maximum)  pressure  mode  shape.  Since  the 
solutions  are  linear,  different  pressure  oscillation  levels  can  be  accounted 
for  by  direct  multiplication  of  the  1.0  psi  results;  e.g.,  for  10  psi, 
multiply  by  10.  The  acoustic  natural  modes  and  frequencies  used  in  the 
present  work  were  supplied  by  Mr.  C.  Bicker  of  NWC. 

The  second  section  of  this  report  gives  a general  overview  of  the 
general  approach  used  in  the  analyses.  The  third  section  contains  a de- 
tailed discussion  of  the  theory  upon  which  the  analysis  is  based  and  presents 
all  applicable  equations.  Details  associated  with  the  NASTRAN  computer 
solution  are  discussed  in  Section  IV.  The  final  two  sections  cover  a dis- 
cussion of  results  and  conclusions.  In  Section  V,  it  should  be  emphasized 
again  that  force  and  displacement  values  given  in  the  text  are  for  a pressure 
oscillation  level  of  1.0  psi  unless  otherwise  stated.  Since  pressure  oscil- 
lation levels  of  10  to  20  psi  or  higher  are  conceivable,  the  values  given  in 
the  text  should  be  multiplied  by  10  or  20  to  obtain  probable  maximum  values. 
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SECTION  II 


GENERAL  APPROACH 


The  general  approach  consisted  of  following  procedures  used  in  a 
previous  program*.  Finite  element  models  were  used  to  represent  the 
entire  space  shuttle  structure.  Mass  and  stiffness  matrices  for  half 
models  of  the  orbiter  and  the  external  tank  with  both  symmetric  and 
antisymmetric  boundary  conditions  were  supplied  by  North  American 
Rockwell  Space  Division  at  Downey,  California.  In  addition,  Rockwell 
furnished  a model,  (mass  and  stiffness  matrices),  for  an  SRB.  A separate 
model  of  the  SRB  structure  above  the  solid  motor  was  provided  by  the 
Marshall  Space  Flight  Center  at  Huntsville,  Alabama.  Hercules  constructed 
a detailed  finite  element  model  of  the  solid  motor  as  a part  of  the  effort 
on  this  program. 

The  NASTRAN  program  on  an  IBM  370/155  computer  was  the  basic  analysis 
tool  used  in  the  program.  Two  versions  of  the  NASTRAN  program  were  used: 

(1)  A NASA  level  15.1  version,  and  (2)  a MacNeal-Schwendler  Company  (MSC) 
program  that  is  approximately  equivalent  to  level  15.5.  The  MSC  version 
contains  a cyclic  symmetry  option  in  the  Frequency  Response  Rigid  Format, 
(R.F.  8).  The  cyclic  symmetry  capability  in  MSC  NASTRAN  was  used  to  analyze 
the  Hercules  solid  rocket  motor  (SRM)  model. 

Attaching  the  models  together  to  obtain  a single  large  model  for  the 
total  shuttle  vehicle  would  have  violated  the  conditions  that  allow  the  use 
of  a cyclic  symmetry  model.  Therefore,  a mechanical  impedance  type  approach 
was  used  which  allowed  each  model  to  be  analyzed  separately.  Results  from 
the  separate  analyses  were  then  combined  to  obtain  the  response  of  the  total 
assembled  space  shuttle  vehicle.  Details  of  the  approach  are  given  in  the 
following  section. 


IF.  R.  Jensen,  Analytical  Prediction  of  Motor  Component  Vibrations 
Driven  by  Acoustic  Combustion  Instability.  Final  Report  AFRPL-TR- 
76-11,  Hercules  Incorporated,  for  the  Air  Force  Rocket  Propulsion 
Laboratory,  Edwards,  CA,  February  1976. 


SECTION  III 


DETAILED  ANALYSIS  METHOD 


The  purpose  of  this  section  of  the  report  is  to  provide  details  on 
how  the  analysis  was  performed.  The  various  finite  element  models  are 
described  in  detail,  general  application  of  mechanical  impedance  is  discussed, 
and  the  aquations  and  other  details  applicable  to  the  shuttle  analysis  are 
provided . 

A.  FREQUENCY  RESPONSE  ANALYSES 

For  a finite  element  model  with  viscous  damping,  the  equations  of 
motion  are: 

[M]  jU|  + [B]jU)  + [K j ju  j = | F ( t ) j (1) 

where:  [M]  = the  mass  matrix 

[B]  = the  viscous  damping  matrix 

[K]  = the  stiffness  matrix 

ju)  = the  displacement  vector 

|F(t) ) = the  applied  load  vector 

For  a harmonic  forcing  function  at  a particular  frequency,  such  as 

jF(t)  ) = jFje*-U)t,  the  equations  representing  the  steady  state  motion  2 are: 

(-ui2[M]  + iuj[B]  + [K])  jU  j = (F ) (2) 

A common  method  of  handling  the  damping  is  to  assume  that  elements  in 
the  damping  matrix  are  proportional  to  corresponding  elements  in  the  stiffness 
matrix3 


[B]  = (g/u.')  [K] 


Equation  (2)  then  becomes 


(V[M]  + (1  + ig)[K])  ju)  = j F } 


2 The  NASTRAN  Theoretical  Manual,  (Level  15) , R.  H.  MacNeal,  Ed.,  April  1972, 
NASA  SP-221 (01) , NASA,  Washington,  D.C.,  page  12.1-3. 

3 Ibid.,  p9.3-8. 


•; 


j 
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Equation  (3>  is  solved  by  the  NASTRAN  program  when  analyses  are  performed 
using  the  frequency  response  rigid  format^.  Let 

[D]  = (V[M]  + (1  + ig)[Kl) 

Then  equation  (3)  can  be  written 

[D]iU|  = {F| 

and  the  solution  to  equation  (3)  can  be  written 

|U|  = [D-l]lF|  (4) 

Let 

[R]  = [D‘l] 

Then  equation  (4)  can  be  written  as 

\U ! = [R]1FJ  (5) 

Using  the  terminology  suggested  in  Reference  5 , [Rj  is  called  a receptance 
matrix . 

Where  a receptance  matrix  is  desired  for  only  a fraction  of  the  degrees 
of  freedom  in  the  finite  element  model,  the  displacement  vector  can  be 
partitioned 

CM  = htiMi  fill 

\ U2J  [^21  I R22J  ( F2  I 

Solving  for  with  ) F 2 1 = 0: 

lui!  = [Rn]!Fii 

This  method  of  obtaining  a reduced  receptance  matrix  is  not  efficient 
because  the  large  D matrix  must  be  inverted  to  obtain  the  R matrix. 


(6) 


(7) 


^The  NASTRAN  User's  Manual,  (Level  15) , C.  W.  McCormick,  Ed.,  June  1972, 
NASA  SP-222(01) , NASA,  Washington,  D.C.,  p3.9-ll. 

^Bishop,  R.  E.  D.,  and  Johnson,  D.  C.,  The  Mechanics  of  Vibration, 
Cambridge  at  the  University  Press,  1960,  London,  England. 
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A more  efficient  way  to  obtain  a reduced  R matrix  is  to  apply  unit 
loads  at  each  of  the  coordinates  in  F3  and  solve  for  Uj,  using  the  regular 
NASTRAN  R.F.-8  solution.  To  clarify  this  approach,  equation  (7)  can  be 
written  as  follows: 


/ > 

f.  ^ 

/ U1 

rll 

r 12 

r 13 

— r In 

fl 

u2 

r21 

r22 

r23 

r2n 

f2 

\ U3 

r 31 

r32 

*33 

— r 3n 

f3 

1 

— 

> 

«n 

V J 

rnl 

rn2 

rn3 

rnn 

or : 

Ul  = rL1fi  + r 12f 2 + r13f3  + ---  + rinfn 
U2  = t2lf l + r22f2  + r23f3  + ---  + r2nfn 
U3  = r3ifl  + r32f2  + r33f 3 + ---  + r3nfn 


Un  rnlfi  + rn2^2  + rn3^3  + + rnn^n 

In  the  above  equations,  when  f^  = 1.0  and  f^  = 0,  i / 1,  the  solution  for 
the  (U3)'s  gives  the  first  column  in  the  desired  reduced  receptance  matrix: 

r 11  = U1 
r2 1 = u2 
r 31  = u3 

rnl  = Un 

The  other  columns  in  the  receptance  matrix  can  be  determined  by  solving 
for  other  force  terms  of  unit  value.  For  example,  a solution  with  f2  = 1.0 
and  all  other  f^  = 0 would  provide  the  second  column  in  the  receptance  matrix. 

A DMAP  ALTER  must  be  used  Co  save  the  displacements  calculated  by 
NASTRAN: 


ALTER  141 

OUTPUT 2 UDVF , , , ,//C,N,-l/C,N,17/C,N,  UDVFTP$ 

END  ALTER 

When  the  above  ALTER  is  used  in  the  executive  control  deck,  the  displacement 
matrix  UDVF  is  written  on  tape  unit  17  and  can  be  saved  for  later  use.  In 
the  same  computer  run  or  in  another  run,  UDVF  can  be  partitioned  to  obtain 
the  R matrix. 

In  the  space  shuttle  analysis,  the  receptance  matrix  at  the  inter- 
connection coordinates  was  required  for  each  model.  The  appropriate 
receptance  matrices  were  obtained  by  applying  unit  forces  at  each  connection 
coordinate  for  each  finite  element  model. 
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B. 


MECHANICAL  IMPEDANCE  TECHNIQUES 


"Impedance"  and  "admittance"  are  terms  generally  associated  with 
electrical  circuits.  The  terms  "mechanical  impedance"  and  "mechanical 
admittance"  are  normally  used  to  indicate  that  an  analogy  is  being  made 
between  an  electrical  circuit  and  a mechanical  system.  The  literature 
on  mechanical  vibration  analysis  contains  a large  amount  of  information  on 
mechanical  impedance-type  approaches.  For  example,  the  Shock  and  Vibra- 
tion Bulletin  contains  many  papers  on  application  of  mechanical  impedance 
techniques . b 

Mechanical  impedance  is  a ratio  of  force  to  velocity.  Mechanical 
admittance,  commonly  called  "mobility,"  is  the  inverse  of  mechanical 
impedance,  i.e.,  a ratio  of  velocity  to  force.  A basic  discussion  on 
mechanical  impedance  and  mobility  can  be  found  in  Reference  7.  The  term 
"receptance"  is  used  to  denote  the  ratio  of  displacement  to  force.  The 
concept  of  receptance  is  discussed  in  References  7,  8,  and  9.  Additional 
discussion  on  electromechanical  analogies  are  contained  in  References  10 
and  11. 

The  term  "immittance"  has  been  used  to  represent  impedance  or  admittance. 
Mechanical  immittance  and  transmission  matrix  concepts  are  discussed  in 
References  12,  13,  and  14. 


b Index  to  the  Shock  and  Vibration  Bulletins.  February  1968,  The  Shock  and 
Vibration  Information  Center,  Naval  Research  Laboratory,  Washington,  D.C. 

^Harris,  C.  M. , and  Crede,  C.  E.,  Shock  and  Vibration  Handbook.  Vol.  1, 
Chapter  10,  McGraw-Hill  Book  Co.,  New  York,  1961. 

O 

Bishop,  R.  E.  D.,  Gladwell,  G.  M.  L.,  and  Michaelson,  S.,  The  Matrix 
Analysis  of  Vibration,  Section  5.5,  Cambridge  at  the  University  Press, 
London,  1965. 

9 

Bishop,  R.  E.  D.,  and  Johnson,  D.  C.,  The  Mechanics  of  Vibration. 
Cambridge  at  the  University  Press,  London,  1960. 

^ Craf ton,  P.  A.,  Shock  and  Vibration  in  Linear  Systems,  Harper  and 
Brothers,  New  York,  1961. 

^MacNeal,  R.  H.,  Electric  Circuit  Analogies  for  Elastic  Structures, 

Vol  2,  John  Wiley  and  Sons,  New  York,  1962. 

12 

Rubin,  S.,  Review  of  Mechanical  Immittance  and  Transmission  Concepts, 
Presented  at  the  71st  Meeting  of  the  Acoustical  Society  of  America, 
Boston,  Mass.,  June  1966. 

13 

JRubin,  S.,  Class  Notes  distributed  at  UCLA  Short  Course  on  Structural 
Dynamics  Analysis,  Los  Angeles,  California,  1967. 

^ Rubin,  S.,  On  the  Use  of  Eight-Pole  Parameters  for  Analysis  of  Beam 
Systems,  Soc . of  Automotive  Engineers,  Reprint  925F , October  1964. 
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When  a sinusoidal  force  drives  a linear  system,  the  steady  state 
response  displacements,  velocities,  and  accelerations  are  sinusoidal  at  the 
frequency  of  the  driving  force.  For  a damped  system,  the  response  is  out- 
of-phase  with  the  driving  force.  The  relationship  between  driving  force 
and  response  can  be  expressed  by  algebraic  equations  involving  complex 
numbers,  (such  as  equation  5).  The  analysis  of  such  a system  is  called  a 
"frequency  response  analysis."  The  use  of  frequency  response -type  analyses 
is  implied  when  mechanical  impedance  is  discussed. 

The  concept  of  impedance  or  receptance  of  the  space  shuttle  structure 
at  a set  of  points  was  used  to  allow  the  individual  shuttle  models  to  be 
analyzed  separately.  The  results  from  the  individual  analyses  were  combined 
to  obtain  solutions  that  represent  the  total  structure. 

Equation  (7)  represents  a linear  system.  Since  the  system  is  linear, 
the  response  for  two  different  load  sets  applied  simultaneously  can  be 
determined  by  applying  each  load  set  separately  and  summing  the  results. 
Consider  one  SRB  undergoing  acoustic  pressure  oscillations.  The  SRB  would 
be  subjected  to  two  separate  loading  systems: 

(1)  The  acoustic  natural  mode  would  load  the  solid  motor  by  means 
of  a certain  pressure  distribution  in  the  motor  combustion 
cavity . 

(2)  The  remainder  of  the  space  shuttle  vehicle  would  apply  loads  to 
the  SRB  at  the  SRB  attach  points  as  the  total  vehicle  vibrates 
in  response  to  the  pressure  oscillations. 

The  objective  of  this  work  is  to  calculate  the  second  load  set. 


In  the  NASTRAN  finite  element  models,  a certain  set  of  displacement 
coordinates  represents  the  SRB  attach  coordinates.  The  SRB  is  attached  to 
the  ET  at  nodes  303,  310,  and  311.  The  node  locations  and  x,  y,  z 
coordinate  directions  are  defined  in  a later  section  of  this  report.  The 
attach  point  displacement  coordinates  are: 


lui 


"u  30  3x 
U 30  3y 
Ju303z  , 
u310y 
u3l0z 
tU  31  ly . 


(8) 


u 

r» 

t : 

£ 


For  this  linear  system,  the  total  response  at  the  attach  points,  |U^j,  can 
be  obtained  by  summing  the  responses  due  to  the  two  separate  load  sets 
discussed  above 

iUTi  = KIsrb  + *Uc!  <9> 


where:  t^oisRB  = t*le  resPonse  at  attach  coordinates  due  only  to  the 

acoustic  pressure  mode 

jUcj  = the  response  at  the  attach  coordinates  due  to  the 

attach  loads  applied  by  the  remaining  shuttle  structure, 
i.e.  the  total  shuttle  less  one  SRB. 

The  displacement  response  ju0j  can  be  calculated  directly  by  using  a 
NASTRAN  model  of  an  SRB.  To  obtain  jUQJ,  a cyclic  symmetry  model  of  an  SRM 
was  analyzed  to  determine  response  to  the  first  acoustic  natural  mode  at 
15.25  Hz.  Details  of  the  jU0)  calculation  are  discussed  in  a following 
section. 

If  the  receptance  of  the  SRB  at  the  attach  points  is  denoted  [Rsrb3 > 
then  | U c ! can  be  expressed  as 

|UC}  = [Rsrb3|fc|  (10) 

where:  jFc)  = the  set  of  forces  applied  to  the  SRB  at  the  SRB  attach 

points 

By  way  of  further  explanation  of  |FCJ,  the  set  of  forces  jFcj  are  internal 
forces  that  occur  at  the  attach  points  between  the  SRB  and  the  ET  due  to 
pressure  oscillations  in  the  SRM.  Cutting  the  structure  at  the  attach  point 
to  show  free  body  diagrams  would  yield  a diagram  showing  internal  forces 
|Fcj  applied  to  the  SRB  and  equal,  but  opposite,  forces,  { -Fc j , applied 
to  the  ET. 

An  equation  similar  to  equation  (10)  can  be  written  for  the  remainder 
of  the  space  shuttle  vehicle.  When  one  SRB  has  been  removed,  the  remaining 
structure  consists  of  the  ET,  the  orbiter,  and  the  other  SRB.  If  the 
receptance  of  the  remaining  structure  at  the  attach  points  is  denoted 
[RrSS],  then  the  forces  applied  result  in  ju^ j displacements 

1ut 1 = [RrssII-Fc!  (11) 

For  this  analysis,  the  only  forces  applied  to  the  remaining  shuttle 
structure  are  the  forces  j-Fcj;  therefore,  the  displacements  in  equation 
(11)  are  the  total  displacements,  jUfj,  rather  than  the  displacements  due 
to  connection  forces,  |UC } > as  in  equation  (10).  Equation  (11)  can  be 
solved  for  connection  forces  jFcj  and  the  result  substituted  into  equation 
(10)  to  eliminate  the  unknown  forces 

jUci  = -[Rsrb][Rrss3IUtI  (12) 

When  equation  (12)  is  substituted  into  equation  (9) , then  the  total 
displacements  are  found  to  be 

juTi  = ([i]  + [RSrb3[Rrss  3) ~ 1 |n0 ! SRB  (13) 
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The  solution  to  equation  (13)  represents  the  desired  response  of  the 
shuttle  to  acoustic  oscillations.  The  connection  forces  can  be  recovered 
by  using  the  solution  of  equation  (11) 

iFc|  = '[rrss^Iut! 

An  equation  similar  to  that  of  equation  (13)  was  used  in  the  component 
vibration  program.  The  Component  Vibration  Final  Report^  contains  a 
simple  example  showing  that  the  technique  of  equation  (13)  is  applicable. 

C.  SPECIFIC  SHUTTLE  ANALYSIS  DETAILS 

To  obtain  the  desired  solution  to  equation  (13)  the  matrices  LRSRB^ > 
rRRss3>  and  iuo!  ggg  must  be  available.  The  calculation  of  these  three 
matrices  is  discussed  below.  The  calculation  of  each  of  the  three  matrices 
is  somewhat  complicated  by  the  fact  that  each  depends  on  solutions  from 
two  different  finite  element  models.  Two  different  models  were  used  to 
represent  the  SRB,  a cyclic  symmetry  model  of  the  SRM  and  a model  of  the 
nose  section  above  the  SRM.  Models  of  the  ET  and  the  orbiter  with  symmetric 
and  antisymmetric  boundary  conditions  are  used  to  represent  the  remaining 
shuttle  structure  for  calculation  of  [Rrss]* 

Locations  of  the  nodes  that  represent  the  attach  points  are  indicated 
in  Figure  1.  The  figure  also  shows  the  basic  X,  Y,  Z coordinate  system 
used  throughout  this  report.  The  coordinate  system  shown  in  Figure  1 was 
adopted  because  it  coincides  with  the  coordinate  system  used  for  definition 
of  the  ET,  orbiter,  and  SRB  models  furnished  by  North  American  Rockwell. 

The  actual  mass  and  stiffness  matrices  involved  were  transmitted  to  Hercules 
on  a computer  tape. 


The  calculation  of  [Rggg]  is  discussed  first.  Figure  2 shows  a 
sketch  of  one  SRB  divided  into  two  parts,  an  SRM  and  the  structure  above 
the  SRM.  For  purposes  of  discussion,  the  structure  above  the  SRM  is 
referred  to  as  the  nose  cone.  This  use  of  "nose  cone"  is  not  in  agreement 
with  official  NASA  nomenclature.  The  SRM  was  represented  by  a cyclic 
symmetry  NASTRAN  model.  Since  10  degree  sections  were  used  in  the  cyclic 
symmetry  model,  there  are  36  nodes  around  the  circumference  of  the  model  as 
shown  in  Figure  3.  The  cyclic  symmetry  SRM  is  described  in  a later  section. 

A NASTRAN  cylindrical  local  coordinate  system,  R,  0,  Z,  was  used  to 
define  node  locations  and  displacement  components  in  the  cyclic  symmetry 
model.  The  Z axis  of  the  local  system  is  parallel  to  the  X axis  of  the 
global  (Rockwell)  coordinate  system.  The  relationship  between  the  local 
R,  0 and  the  Global  Y,  Z axes  is  indicated  in  Figure  3.  Attach  nodes  310 
and  311  are  located  circumferentially  as  shown  in  Figure  3.  The  transforma- 
tion below  was  written  to  relate  the  displacements  in  the  two  coordinate  systems. 


^ Analytical  Prediction  of  Motor  Component  Vibrations  Driven  by  Acoustic 
Combustion  Instability,  op.  cit. 
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NODE  302. 


NODE  304. 


NODE  301 


FORWARD  SRB  ATTACH  POINT,  NODE  303 


NOSE  CONE  (SEE  TEXT) 


^Ui^NC’  (Fi^NC 

Ui,  Fi 


FORWARD  SKIRT 


AFT  SRB  ATTACH  POINTS,  NODES  310  & 311 


Figure  2.  Sketch  of  the  SRB  Showing  the  Division  into  SRM 
and  Nose  Cone  Models 
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Using  compact  notation 


iU310,  nl  “ CT]|U310)  uj  (15) 

Also 

lU310,  111  = [T1]  1u310 > u!  (15a) 

The  unprimed  displacement  vector  in  the  Y and  Z directions  contains  only 
three  components  because  the  SRB  is  attached  to  the  ET  in  only  those  three 
components  at  nodes  310  and  311. 

The  cyclic  symmetry  model  was  used  to  obtain  a receptance  matrix  for 
the  SRM  by  applying  unit  loads  at  all  connection  points,  as  explained 
previously.  The  nodes  numbered  1 through  36  in  Figure  3 are  located  at 
the  forward  end  of  the  forward  skirt  of  the  SRM  and  are  therefore  the 
connection  points  between  the  SRM  and  the  nose  cone.  The  receptance  matrix 
obtained  from  analysis  of  the  cyclic  symmetry  model  had  the  form  shown  below. 


The  receptance  matrix  from  the  SRM  analysis  was  partitioned  in 
equation  (16)  to  separate  the  L'i  displacements  from  the  attach  point 
displacements.  Equation  (16)  can  be  expanded  as  follows 

Ui  = RsH  Fi  + rS12  f310,  11  (17a) 

U310,  11  = RS21  Fi  + RS22  F310,  11  (17b) 

The  transformation  of  equation  (15)  also  holds  when  forces  are  substituted 
for  displacements 

If310,  111  * [T]!F3io,  ni  (18> 

Substituting  equations  (15)  and  (18)  into  equations  (17a)  and  (17b)  gives 

= RS11  Fi  + (RS12  T>F3io,  11  (19a) 

T U310,  11  = RS21  Fi  + RS22  ^ F310,  11  (19b) 

tr  

Since  multiplying  (T  ) by  (T)  gives  the  identity  matrix,  equation  (19b)  can 
be  rewritten  as 

u310,  11  = (^  RS2l)pi  + (^  RS22  F)f310,  11  (19c> 


The  nose  cone  model  only  contains  four  nodes  as  shown  in  Figure  2. 
Details  on  how  the  nose  cone  model  was  obtained  are  given  in  a later  section. 
Unit  forces  were  applied  at  nodes  303  and  304  to  obtain  the  nose  cone 


r 


Again,  the  receptance  matrix  is  partitioned  to  separate  the  ET/SRB  attach 
degrees  of  freedom  from  the  SRM/nose  cone  interface  degrees  of  freedom. 
Equation  (20)  can  be  expanded  as  follows: 


303 

RN11  F303  + RN12  FiNC 

(21a) 

iNC 

= rN21  F303  + RN22  FiNC 

(21b) 

At  this  point,  another  transformation  equation  is  needed  to  relate  the 
displacements  at  36  nodes  on  the  SRM  to  the  U^fjC  displacements  at  the 
single  node  (304)  on  the  nose  cone.  The  angle  Cpi  is  measured  counterclockwise 
from  the  R axis  as  shown  in  Figure  3.  The  required  transformation  for  the 
itb  node  of  the  36  SRM  interconnection  nodes  is  then 

( \ r i (*  \ 

UiR  0 -Sin^+S)  Cos(¥i+5)  0 0 0 u304x  \ 

UieV  = 0 -Cos(cPj+5)  -Sin^+S)  r 0 0 U304Y 

Uiz  1 0 0 0 rCos (<Pj+5)  rSin(cpt+5)  <U304z  l 

^ ' L J 

U 304RX 
U304RY 
^U304RZ 

When  the  above  equation  is  written  to  include  all  36  nodes,  the  following 
compact  expression  is  used 

l»i!  = [T]jUiNC!  (22) 

The  internal  forces  applied  to  the  SRM  are  equal  in  magnitude  but 
opposite  in  direction  from  those  applied  to  the  Nose  Cone.  Therefore,  the 
force  transformation  corresponding  to  equation  (22)  is: 

tFf!  = -[T]|FiNC!  (23) 

Substituting  equations  (22)  and  (23)  into  equations  (19a)  and  (19c)  gives 

T UiNC  = ‘RSll  T FiNC  + RS12  * F310,  11  (24a> 

U310,  11  = 'TT  rs21  T FiNC  + ^ RS22  F F310,  11  (24b> 

Substituting  equation  (21b)  into  equation  (24a)  and  rearranging  gives 

(T  RN22  + Rsll  T^FiNC  = rS12  * F310,  11  ' T RN21  F303 
Premultiply  by  the  transpose  of  T 

TT(T  RN22  + Rsll  T^FiNC  = tT  RSl2  ? F310,  11  * tT  t rN21  F303 

Let 

a = TT(T  Rn22  + Rsn  T)  (25) 
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Then 


FiNC  = a'1  Tt  RS12  T F310j  n - a"1  TT  T RN21  F303  (26) 

Equation  (26)  can  be  used  to  eliminate  Fj^q  from  equations  (21a)  and  (24b) 
U303  = (Full  ‘ rn120'  Li'TtRN21^F303  + RN120'  1tTr S 1 2TF 3 10 , 11  (27a) 

U310,  11  = TTRS21Ta‘lTlraN21F303  + ^22?  * TTrS21T“'1tTrs12T>F310,  u 

(27b) 

Equations  (27a)  and  (27b)  can  be  combined  to  give  an  expression  for  the  SRB 
receptance  matrix 

F303 


F 310 , ll| 
(28) 

The  square  matrix  in  equation  (28)  is  the  receptance  matrix  for  one  SRB, 

LRSRBJ  • 

The  calculation  of  |l'0}gRg  is  discussed  below.  JUo!gRg  is  the  response 
of  the  SRB  at  the  attach  points  to  a particular  acoustic  natural  mode.  The 
calculation  of  {U0|SRB  similar  to  the  calculation  of  [Rsrm]  discussed 
above.  The  response  cannot  be  determined  by  direct  analysis  because  one 
attach  point  is  on  the  nose  cone  while  the  other  two  attach  points  are  on 
the  SRM. 

Considering  the  SRM  to  be  isolated  from  the  nose  cone,  the  response 
at  the  SRM  attach  points,  nodes  310  and  311,  is  obtained  by  adding  the 
response  due  to  acoustic  mode  (p)  forces  to  the  response  due  to  internal 
(Ft)  forces  at  the  SRM/nose  cone  interface 

iu3io,  h!t  = lu3io , nip  + iu3io,  11  If  (29) 

The  response  at  the  SRM/nose  cone  interface  is  expressed  in  a similar 
manner : 


\ U303 

L 

(%11  - RN12Q  T TRN21>  I (RNUa'^RgizT) 

1 

[_U310,  11 

r 

_(TTRs21Tcr1TTTRN21)  j (TTRS22T  - TTRs21Tcv'1TTRsl2T) 

|UiiT  = |UiSp  + lUilp  (30) 

The  response  of  the  SRM  to  the  acoustic  mode  is  obtained  directly  by 
applying  loads  that  represent  the  acoustic  mode  pressure  distribution  to  the 
NASTRAN  cyclic  symmetry  model.  Therefore,  jU3xo,  11*D  anc*  !ui!n  are  obtained 
by  partitioning  the  UDVF  displacement  matrix  from  a NASTRAN  solution.  Note 
that  equations  (29)  and  (30)  are  similar  to  equation  (9). 
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The  terms  jUg^o  11 !f  and  I ^ i I F represent  the  response  of  the  SRM 
to  the  interface  forces,  Fi.  Using  equations  (19a)  and  (19c)  with 
F310  33  = 0,  the  following  expressions  are  obtained 


(Ui)p  = Rgi3  F 3 (31a) 

(U310,  11>F  = (TT  RS2l)Fi  (31b) 

Substituting  (31a)  and  (31b)  into  equations  (29)  and  (30)  gives: 

<U310,  11>T  = 0*310,  ll)p  + (T1  RS2l)Fi  (32a) 

(Ui)T  = (Ui)p  + Rsn  F i (32b) 

Substituting  equation  (23)  into  (32a)  and  substituting  equations  (22)  and 
(23)  into  equation  (32b)  results  in 

(u310,  lpT  = (u310,  lPp  * ^ RS21^T  FiNC  (33a) 

T(UiNC)x  = (U£)p  - Rsn  T FiNC  (33b) 

For  the  nose  cone,  only  the  response  to  interface  forces  iFiNc!  is  required . 
Using  equations  (21a)  and  (21b)  with  F303  equal  to  zero  gives 

(0303)7  = rN12  FiNC  (34a) 

(UiNC^T  = rN22  FiNC  (34b) 


Substituting  equations  (34b)  into  equation  (33b)  premultiplying  by  TT  and 
rearranging  gives 

TT(RS11  T + TRN22)FiNc  = TX(Ui)p 

The  coefficient  of  Fj^c  in  the  last  equation  can  be  recognized  as  a,  as 
defined  in  equation  (25). 

Therefore , 

FiNC  = TT«i>P  (35) 

Substituting  equation  (35)  into  equations  (33a)  and  (34a)  gives 

0*310,  11>T  = (u310,  ll>p  - (TTrs21>T  (36a> 

(u303)t  “ Rsi^’^OJiJp  (36b> 

In  equation  (36a) , equation  (15a)  is  substituted  to  obtain  the  result  in 
the  primed  system  that  corresponds  to  direct  NASTRAN  output: 

(1*310 , 11>T  = ^(0310,  ll>p  - (TTRS2l)T  «'LrT(Ui)p  (36c) 
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Using  equations  (36b)  and  (36c),  an  expression  can  be  written  for  jU0jsRg 

(U 30 3> x f(RN12  a'1  TT)  J 0 f(Ui)p 

1uo!srb  = ~'iz I r (3?) 

((W310,  n)TJ  L(‘^Rs2iT  “ T > : T JI/U310,  ii)p 

The  remainder  of  the  discussion  in  this  section  is  concerned  with  the 
calculation  of  [Rggsl  • Equation  (11)  defines  the  displacement  and  force 
vectors  that  are  associated  with  receptance  matrix  [Rggs7,  The  jU.j.|  vector 
in  equation  (11)  has  individual  components  as  defined  in  equation  (8) . 

NASTRAN  finite  element  models  of  the  ET  and  the  orbiter  were  supplied 
by  the  North  American  Rockwell  Company,  Space  Division,  at  Downey,  California. 
The  furnished  models  represent  only  one-half  of  the  structure  and  separate 
models  were  supplied  for  symmetric  and  asymmetric  boundary  conditions.  The 
X-Z  plane  as  shown  in  Figure  1 was  taken  as  the  plane  of  symmetry  for  the 
models.  When  the  models  with  symmetry  boundary  conditions  are  analyzed, 
the  results  should  represent  a symmetric  structure  subjected  to  symmetric 
loads.  Therefore,  a solution  for  a particular  acoustic  mode  using  the  models 
with  symmetric  boundary  conditions  should  represent  the  condition  where 
both  solid  rocket  motors  are  being  subjected  to  unstable  acoustic  oscilla- 
tions that  are  in-phase.  Use  of  the  models  with  asymmetric  boundary 
conditions  would  represent  the  condition  where  both  SRM' s were  oscillating 
out-of-phase  with  one  another.  The  difference  between  the  symmetric 
and  asymmetric  solutions  should  give  results  for  the  situation  where  only 
one  SRM  is  undergoing  unstable  acoustic  oscillations. 


The  receptance  equation  for  the  ET  is  as  follows 

’u  [rh  J Riil/* F ) 

jUETj  = ---l  = — H — - \-l-\  (38) 

u J r2i  j R22J  |^F  J 

The  juj  vector  in  equation  (38)  has  components  as  defined  in  equation  (8) 
for  both  symmetric  and  asymmetric  solutions.  The  jui  vector  of  equation  (38) 
has  different  components  for  the  symmetric  and  asymmetric  models  as  follows 


SECTION  IV 


NASTRAN  COMPUTER  SOLUTION 


The  preceding  section  (Section  III)  contains  a discussion  of  the 
theory  upon  which  this  analysis  is  based.  In  this  sect. on,  application  of 
the  theory  is  discussed.  The  NASTRAN  computer  program  was  used  as  the  basic 
analysis  tool.  The  cyclic  symmetry  analysis  of  the  SRM  model  and  the  DMAP 
program  analysis  of  the  Rockwell  furnished  models  are  furnished  under 
separate  headings. 

SRM  Analysis  Using  the  Cyclic  Symmetry  Model 

Recent  versions  of  NASTRAN  (i.e.  Level  15.5  and  Level  16.0)  contain  the 
capability  to  analyze  cyclic  symmetric  structures  using  the  static  analysis 
(R.F.-l)  or  real  eigenvalue  analysis  (R.F.-3)  rigid  iormats.  Corresponding 
NASTRAN  documentation16  contains  a description  of  cyclic  symmetry;  the  des- 
cription will  not  be  repeated  here. 

A special  version  of  NASTRAN  was  use'’  for  the  cyclic  symmetry  analysis 
of  the  SRM  model.  The  Mac  Neal-Schwendl •_  ’ Company  (MSC) , working  under  contract 
with  Hercules  Incorporated,  added  the  cyclic  symmetry  capability  to  the  fre- 
quency response  rigid  format,  (R.F.-8),  in  NASTRAN.  The  work  was  sponsored  by 
the  Air  Force  Rocket  Propulsion  Laboratory  at  Edwards  AFB  under  contract 
F0461I-73-C-0025  with  Hercules.  The  MSC  version  of  NASTRAN  can  be  obtained 
from  the  Mac  Neal-Schwendler  Corporation  at  7442  No.  Figueroa  Street,  Los 
Angeles,  California  (90041),  or  it  can  be  used  on  some  of  the  large  computer 
systems  that  lease  it  from  MSC;  e.g.,  the  CD C Cybernet  system. 

The  cyclic  symmetry  analysis  capability  in  NASTRAN  allows  an  efficient 
general  three-dimensional  analysis  to  bi  performed  on  a structure  that  is 
cyclic  symmetric  by  modeling  only  a portion  of  the  structure.  For  a rocket 
motor  with  the  usual  slotted  grain  design,  a radial-axial  plane  passed  through 
the  center  of  each  slot  divides  the  motor  into  sections  which  repeat  around 
the  circumference  of  the  motor.  Since  each  such  section  is  also  symmetric 
about  a radial-axial  plane  that  would  bisect  it,  NASTRAN  requires  a model 
of  only  one-half  of  a section.  The  SRM  has  11  slots  in  the  grain  design.  A 
half  section  model  would  therefore  cover  1/22  of  the  circumference  or  about 
16.4  degrees.  Previous  work1 7 has  indicated  that  a grid  slice  somewhat  less 
than  15  degrees  is  desirable  because  of  the  better  approximation  of  curved 
surfaces  that  a narrow  slice  would  offer.  To  obtain  a slice  narrower  than 
16.4  degrees  two  different  approaches  could  be  used:  1.)  The  16.4  degree 
motor  section  could  be  modeled  with  a grid  containing  two  8.2  degree  slices, 
or  2.)  a model  could  be  constructed  for  a motor  with  a greater  number  of  slots. 
Option  1.)  results  in  a grid  with  more  degrees  of  freedom  and  a stiffness 
matrix  with  a larger  bandwidth.  A significantly  larger  run  time  would  result 
from  the  use  of  option  1.). 


lfjNASTRAN  Theoretical  Manual,  op.  cit. 

1 7 "Analytical  Prediction  of  Motor  Component  Vibrations...",  op.  cit. 


In  an  effort  to  maintain  reasonable  run  times  on  tbe  computer,  option  2,) 
was  adopted.  A motor  model  with  18  slots  was  used  to  represent  the  actual  motor 
which  has  11  slots.  The  slot  shape  in  a radial-axial  plane  was  not  changed  and 
slot  width  was  adjusted  to  keep  the  volume  of  propellant  in  the  slotted  area 
of  the  motor  approximately  constant  between  the  18  slot  and  11  slot  configurations. 
The  18  slot  model  should  provide  a good  representation  of  the  11  slot  motor  in 
the  most  important  longitudinal  fundamental  mode  since  the  major  structural  effect 
the  slots  have  is  to  reduce  the  hoop  stiffness  of  the  propellant  grain  in  the 
slotted  area.  Hercules  has  obtained  good  results  from  a 2-D  axisymmetric  finite 
element  computer  program  that  utilizes  slot  approximation  elements  that  have 
no  hoop  stiffness. 

The  grid  layout  in  a radial-axial  plane  is  shown  in  Figure  4.  A projected 
view  of  the  grid  is  shown  in  Figure  5.  In  both  Figures  4 and  5,  the  actual  grid 
has  been  rotated  180  degrees  about  the  motor  axis  and  plotted  at  both  the  0 and 
180  degree  locations  to  show  the  motor  outline.  The  grid  is  quite  coarse  com- 
pared with  the  usual  solid  rocket  motor  finite-element  grid.  The  coarse  grid 
is  a result  of  efforts  to  model  a very  large  structure  with  an  adequate  model 
and  still  maintain  reasonable  computer  run  times.  Based  on  past  experience, 
the  grid  should  provide  adequate  results  for  low  frequency  analyses  but  would 
have  questionable  accuracy  at  higher  frequencies.  More  discussion  on  grid 
refinement  is  contained  in  the  Component  Vibration  Program  Final  Report^-®. 

To  avoid  a lengthy  description  of  the  model  geometry,  material  properties, 
etc.,  a copy  of  the  NASTRAN  bulk  data  deck  is  included  as  Table  I.  The  corres- 
ponding case  control  deck  is  shown  (abridged)  in  Table  II.  Only  the  first  and 
last  pages  of  the  case  control  deck  are  shown  in  Table  II  because  intermediate 
subcase  cards  are  very  repetitious.  The  corresponding  executive  control  deck 
listing  is  shown  in  Table  III.  The  alter  statements,  with  the  exception  of 
the  0UTPUT2  statement,  were  required  to  use  the  cyclic  symmetry  option  in  MSC 
NASTRAN.  The  0UTPUT2  alter  was  used  to  write  the  displacement  solutions  on 
a computer  tape  that  could  be  saved  for  later  use. 

The  NASTRAN  solution  represented  by  the  data  of  Tables  I,  II,  and  III  was 
run  to  obtain  data  for  the  SRM  receptance  matrix.  A total  of  10  unit  loads 
were  applied  to  the  SRM  model  to  determine  the  receptance  matrix.  Four  unit 
loads  were  applied  at  the  attach  points,  (load  set  Fj^q  n)>  and  six  unit 
loads  were  applied  at  the  SRM/Nose  Cone  interconnection  points.  Unit  loads 
were  applied  at  points  1 and  2 as  shown  in  Fig-  3.  Points  1 and  2 repre- 
sent the  only  types  of  unique  points  among  the  36  different  interface  nodes 
of  Fig.  3.  Node  1 represents  nodes  that  are  half  way  between  slots  and 

node  2 represents  nodes  that  are  in  line  with  slot  tips.  The  computer  solution 
represented  by  Tables  I,  II,  and  III  was  run  on  the  IBM  370/155  with  1000K  of 
core.  The  CPU  time  for  the  run  was  7.07  hours  and  the  total  run  time  (CPU 
plus  wait)  was  8.07  hours.  The  problem  had  1832  degrees  of  freedom  in  the 
analysis  set.  The  1/36  section  grid  therefore  represents  an  equivalent  full 
360  degree  equivalent  model  with  18  x 1832  = 32,976  degrees  of  freedom. 

Results  from  the  analysis  were  the  displacements  at  all  nodes  for  each  of  the 
10  different  loads.  All  displacements  are  in  data  block  UDVF  which  was  written 
on  an  0UTPUT2  tape. 


18  Op.  Cit. 
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A relatively  long  computer  run  time  was  anticipated  for  the  unit  load 
solutions.  Therefore,  an  effort  was  made  to  reduce  the  anticipated  run  time. 

The  method  investigated  for  a reduction  in  run  time  was  the  use  of  a cyclic 
index,  K,  less  than  K^x.  When  the  Mac-Neal-Schwendler  Company  presented  the 
modified  NASTRAN  program  to  Hercules,  they  stated  that  some  problems  could 
probably  be  solved  with  sufficient  accuracy  with  a cyclic  index  less  than 
the  maximum  value.  During  the  Component  Vibration  Program,  Hercules  investi- 
gated the  possibility  of  reducing  the  cyclic  index  by  running  real  eigenvalue 
solutions.  The  conclusion  was  that  all  K indices  were  required  for  the  parti- 
cular motor  being  studied  to  obtain  all  natural  frequencies  up  to  at  least 
500  Hz  and  probably  for  higher  frequencies  as  well.  Based  on  that  study,  all 
cyclic  indices  were  used  for  the  analyses  performed  during  the  component 
vibration  program. 

For  frequency  response  type  solutions,  a more  applicable  way  to  study 
the  need  for  high  values  of  the  K index  is  to  conduct  frequency  response  type 
solutions  with  various  K values.  One  section  of  the  SRM  model  consisting  of 
15  propellant  elements  and  9 case  elements  was  analyzed  for  K = 3,  K = 6,  and 
K = 9.  The  value  K = 9 is  the  maximum  value  for  the  36  slice  SRM  finite 
element  grid.  Some  results  from  these  three  analyses  are  shown  in  Table  IV. 
Based  on  the  results  of  the  comparison  analyses,  the  decision  was  made  to  use 
a K = 6 in  the  main  analysis.  The  K = 6 results  were  judged  to  be  sufficiently 
accurate  representations  of  the  K = 9 solutions.  Using  fewer  than  the  maximum 
number  of  cyclic  indices  in  a solution  is  similar  to  omitting  some  of  the 
higher  frequency  modes  to  reduce  degrees  of  freedom  of  the  system  when  modal 
coordinates  are  used  to  solve  a problem. 

As  mentioned  above,  only  10  unit  loads  were  applied  for  the  SRM  re- 
ceptance  matrix,  [RsRM]>  calculation  run.  To  obtain  the  RsRM  matrix  more 
directly,  three  loads  would  have  been  applied,  one  at  a time,  to  each  of  the 
36  nodes  at  the  SRM/Nose  cone  interface.  Instead  of  obtaining  3 x 36  = 108 
solutions,  the  results  from  six  unit  loads  at  two  nodes  were  used  to  obtain 
all  required  information.  The  results  of  applying  loads  at  node  1 (Fig.  3) 
could  be  rotated  to  apply  to  loads  at  any  other  odd  numbered  node.  For  example, 
the  radial  response  at  node  six  due  to  a unit  axial  load  at  node  three  would 
be  the  same  as  the  radial  response  at  node  4 with  the  unit  axial  load  applied 
at  node  1.  The  responses  due  to  loads  applied  at  nodes  1 and  2 were  rotated 
as  required  by  using  PARTITION  and  MERGE  operations  in  a DMAP  program.  The 
UDVF  data  block  was  operated  on  in  the  DMAP  program  to  obtain  the  desired 
receptance  matrix.  The  DMAP  program  used  consisted  of  about  300  DMAP  state- 
ments and  the  program  required  three  separate  runs  to  complete  due  to  NASTRAN 
program  limitations.  The  matrix  [Rsrm]  is  a complex  matrix  of  order  112.  The 
large  [RsRm3  matrix  is  not  included  in  this  report. 

A separate  computer  run  was  made  to  calculate  the  {Uo]SRM  vector.  Since 
applied  loads  were  symmetric,  the  run  was  made  with  cyclic  index  K = 0.  Using 
a 500K  core,  the  CPU  time  was  70  minutes  and  the  total  run  time  was  about 
83  minutes.  As  in  the  previous  case,  a EMAP  alter  was  used  to  write  the  UDVF 
data  block  on  tape.  A separate  DMAP  program  was  used  to  form  {Uo]sg^  from  the 


UDVF  data  block.  The  DMAP  program  listing  is  given  in  Table  V.  The  calcu- 
lated {UolsRM  vector  is  shown  in  Table  VI.  Ttie  complex  displacements  shown 
in  Table  VI  represent  the  response  of  the  SRM  cyclic  symmetry  model  to  the 
15.25  Hz  first  longitudinal  acoustic  mode,  (see  equations  31a  and  31b). 

Calculation  of  [Rsrb]  and  {U0}sRB 


The  equations  for  the  calculation  of  [RsrrI  and  {Uol.SRR  are  given  in 
Section  III  of  this  report,  (refer  to  equations  28  and  37).  The  DMAP  program 
used  for  these  calculations  is  shown  in  Table  VII.  The  resulting  data  are 
given  in  Table  VIII. 

During  the  time  when  computer  runs  were  being  made  to  set  up  the  calcu- 
lation of  [Rsrb])  the  series  multiply  and  add  module  (SMPYAD)  was  found  to 
contain  an  error.  The  SMPYAD  module  does  not  work  correctly  with  complex 
matrices.  Therefore,  the  DMAP  programs  all  use  only  the  MPYAD  multiply  and 
add  module. 

The  Nose  Cone  receptance  matrix  was  required  in  the  calculation  of 
[RSrb].  The  Nose  Cone  receptance  matrix  was  obtained  from  the  64  degree-of- 
freedom  SRB  model  furnished  by  Rockwell.  The  DMAP  program  used  to  calculate 
the  Nose  Cone  receptance  matrix  is  shown  in  Table  IX.  The  Nose  Cone  receptance 
matrix,  [Rncone]  ts  presented  in  Table  X. 

Analysis  t the  Rockwell  El  and  Orbitcr  Models 


The  Rockwell  finite  element  models  are  characterized  by  mass  and 
stiffness  matrices.  The  mass  and  stiffness  matrices  for  each  model  were 
transmitted  to  Hercules  on  a computer  tape.'  A FORTRAN  program  was  written 
to  read  the  Rockwell  tape  and  write  a NASTRAN  compatible  tape. 


The  NASTRAN  tape  was  used  as  input  to  the  DMAP  sequence  shown  in  Table  XI. 
The  DMAP  program  calculates  receptance  matrices  by  applying  unit  loads  to  the 
applicable  attachment  coordinates.  The  equations  of  motion  as  defined  in 
equation  3 are  solved  to  determine  displacement  response  to  the  unit  loads. 

A structural  damping  factor  of  g = 0.06  has  been  used  in  the  calculations. 

The  receptance  matrices  obtained  from  the  EMAP  program  are  given  in  Table  XII. 

Calculation  of  the  Force  and  Displacement  Response 

The  displacements  and  forces  at  the  attach  points  were  calculated  bv 
the  DMAP  program  shown  in  Table  XIII,  (Refer  to  equations  13,  14,  42,  and  43b). 
The  output  from  the  program  shown  in  Table  XIII,  is  given  in  Table  XIV.  The 
following  terminology  is  used  in  Table  XIV: 
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Matrix  Name  in  Table 


UTOT 

FTOT 

UBAR 

FBAR 

RRSS 

UTOA 

FTOA 

UBAA 

FBAA 

RRSA 

USIG 

FSIG 

UBSG 

FBSG 


Matrix 

(uT] 

PT) 

(U} 

{F} 

[rrss3 

{ut1 

[ft] 

[v] 

£f] 

[rrss1 

^■UTj1SYM  + ^ UT-  ASYM 

^ft*sym  + ^ft''asym 

^SYM  + ASYM 
^SYM  + '-F^ASYM 


Boundary  Conditions 
Symmetric 
Symmetric 
Symmetric 
Symmetric 
Symmetric 
Antisymmetric 
Antisymmetric 
Antisymmetric 
Antisymmetric 
Antisymmetric 
Not  Applicable 
Not  Applicable 
Not  Applicable 
Not  Applicable 
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SECTION  V 


DISCUSSION  OF  RESULTS 


Response  of  the  Space  Shuttle  Vehicle  to  the  first  longitudinal 
acoustic  mode  has  been  calculated.  A sketch  of  the  vehicle,  indicating 
nodes  at  the  interconnection  points,  is  shown  in  Figure  1.  As  shown  in 
Figure  1,  Node  303  represents  tlie  forward  SRB/F.T  attach  point.  Nodes  310 
and  311  are  the  aft  SRB/ET  attach  points.  Nodes  90  and  91  represent  the 
connections  between  the  ET  and  the  Orbiter.  The  X,  Y,  2 coordinate  system 
shown  in  Figure  I is  the  system  used  with  the  Rockwell/NASA  models  of  the 
ET,  SRB,  and  Orbiter. 

Based  on  symmetry  about  the  X-Z  plane,  the  problem  was  solved  with 
both  symmetri  and  anti-symmetric  boundary  conditions.  The  solution  based 
on  symmetry  boundary  conditions  represents  the  situation  where  both  SRM's 
are  under  going  unstable,  in-phase  pressure  oscillations.  The  solution  based 
on  the  anti-symmetry  boundary  conditions  represents  the  corresponding  out-of- 
phase  situation.  The  sum  of  the  symmetric  and  asymmetric  solutions  represents 
the  situation  where  only  one  SRM  undergoes  unstable  pressure  oscillations. 

The  forces  and  displacements  from  the  in-phase  analysis  are  given  in 
Tables  XV  and  XVI.  The  largest  force  for  the  + I psi  pressure  oscillation 
level  is  1,610  pounds  at  Node  303  in  the  Y direction.  For  a pressure 
oscillation  level  of  + 10  psi,  the  maximum  force  would  increase  to  16,100 
pounds.  The  corresponding  displacements  and  forces  from  the  out -of-phase 
analysis  are  shown  in  Tables  XVII  and  XVIII.  The  maximum  force  from  the 
"out-of-phase"  solution  is  1,012  pounds  at  Node  30 3 in  the  X direction. 

When  the  in-phase  and  out-of-phase  solutions  are  added  to  obtain  the 
results  for  the  situation  where  only  one  booster  is  undergoing  unstable 
pressure  oscillations,  a maximum  interface  force  of  2,127  pounds  is  obtained 
at  Node  303  in  the  X direction.  However,  such  addition  of  the  solutions 
represents  a pressure  oscillation  level  of  + 2 psi.  The  normalized  maximum 
force  for  a + 1 psi  pressure  level  would,  therefore,  be  1,064  pounds.  The 
"in-phase"  condition,  therefore,  produces  the  greatest  forces. 

Although  reasonable  care  was  taken  during  construction  and  checkout  of 
the  finite  element  models  used  in  this  analysis,  and  in  the  set-up  of  the 
analysis  procedure,  the  results  could  be  in  error.  Because  of  the  relatively 
involved  procedure  that  was  used  to  calculate  attach  point  forces,  an  error 
could  be  difficult  to  detect.  Some  kind  of  a check  on  results  would  be 
desirable.  Since  time  and  budget  would  not  allow  for  analysis  of  a very 
simple  model  by  the  detailed  procedure  for  checkout  purposes,  an  attempt 
has  been  made  to  obtain  at  least  a rough  order-of-magnitude  checkout  by 
performing  some  simple  hand  calculations. 
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TABLE  XV 

Displacements  at  the  Attach  Points 
Due  to  Symmetric  Oscillation  in  the  First 
Acoustic  Mode  at  15.25  Hz.  (i  1 psi  pressure  oscillation  level) 


DISPLACEMENT 

COMPONENT 


DISPLACEMENT 
AMPLITUDE  (IN. 


DISPLACEMENT 
PHASE  (DEG 


U303X 

.00094 

96 

U303Y 

.00063 

38 

U303Z 

.00005 

-179 

U310Y 

.00054 

87 

U310Z 

.00024 

97 

U311Y 

.00056 

88 

U90X 

.00039 

-32 

U9OZ 

.00012 

172 

U91X 

.00019 

-6 

U9iY 

.00001 

-117 

4 

1 U91Z 

.00011 

144 

TABLE  XVI 


Forces  at  the  Attach  Points  Due  to 
Symmetric  Oscillation  in  the  First  Acoustic 
Mode  at  15.25  Hz.  (±  1 psi  pressure  oscillation  level) 


FORCE  FORCE  FORCE 

COMPONENT  AMPLITUDE  (LB.)  PHASE  (DEG,). 


F303X 

1306 

95 

F303Y 

1610 

49 

F3032 

28 

174 

f310Y 

391 

91 

f310Z 

89 

-64 

F3L1Y 

356 

93 

F90X 

1 

-32 

F 

90Z 

14 

148 

F91X 

152 

-67 

F9iy 

38 

82 

F91Z 

59 

-104 

72 


IMMIMkfll 


- 


TABLE  XVII 


Displacements  at  the  Attach  Points  Due 
to  Asymmetric  (out-of -phase)  Oscillation  in 
the  First  Acoustic  Mode  at  15.25  Hz. 

(±  1 psi  pressure  oscillation  level) 


DISPLACEMENT 

COMPONENT 

DISPLACEMENT 
AMPLITUDE  (IN.) 

DISPLACEMENT 
PHASE  (DEC.) 

U303X 

.00048 

136 

U303Y 

.00029 

-20 

U303Z 

.00012 

117 

U310Y 

.00044 

95 

U310Z 

.00008 

-98 

U311Y 

.00056 

99 

U90Y 

.00045 

-159 

U91X 

.00001 

139 

U91Y 

.00012 

29 

U91Z 

.00017 

-82 

TABLE  XVIII 


FORCE 

COMPONENT 

F303X 

F303Y 

f303Z 

f310Y 

f310Z 

F311Y 

F90Y 

F9IX 

F91Y 

F91Z 


Forces  at  the  Attach  Points  Due  to 
Asymmetric  (out-of-phase)  Oscillation  in 
the  First  Acoustic  Mode  at  15.25  Hz. 

(±  1 psi  pressure  oscillation  level) 


FORCE  FORCE 

AMPLITUDE  (LB.)  PHASE  (DEG.) 


1012 

142 

249 

51 

64 

117 

507 

-108 

219 

-77 

599 

-108 

150 

93 

184 

81 

235 

89 

83 

85 
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Estimation  of  Longitudinal  Mode  Attach  Forces  by  Hand  Calculation 


To  make  a crude  estimate  of  the  attach  forces,  the  SRB  was  considered 
to  be  a rigid  body.  When  forces  representing  the  chamber  pressure  are 
applied  to  the  SRM  finite  element  model,  a net  axial  force  on  the  model 
is  produced.  The  net  force  occurs  because  the  nozzle  opening  in  the  aft 
dome  does  not  carry  any  load  to  offset  the  force  produced  on  a correspond- 
ing area  in  the  forward  dome.  If  the  pressure  applied  to  the  finite  ele- 
ment model  were  an  oscillatory  pressure,  then  the  net  force  would  be 
oscillatory.  The  net  force  thus  applied  would  cause  oscillatory  axial 
accelerations  of  the  total  SRB  considered  as  a rigid  body. 

A uniform  internal  static  pressure  was  applied  to  the  SRM  finite  ele- 
ment model.  The  model  was  constrained  in  the  axial  direction  at  the  SRM/ 
Nose  Cone  attach  points.  For  the  10  degree  slice  model  with  symmetry 
boundary  conditions  applied  at  the  slice  sides,  (i.e.,  radial-axial  planes 
axial  constraint  forces  of  49.58  lbs  were  calculated  for  each  side  of  the 
grid,  a total  of  99.16  lbs  for  the  10  degree  slice.  For  a complete  360 
degree  motor  the  1 psi  thrust  would  therefore  be  36  x 99.16  = 3569.76  lbs. 
For  a chamber  pressure  of  850  psi  the  thrust  of  850  x 3569.76  = 3.03  x 
10^  lbs  compares  well  with  the  3.1  x 10^  lbs  estimated  by  Space  Shuttle 
engineers.  This  good  agreement  gives  us  confidence  that  the  SRM  model 
is  yielding  reasonable  rigid  body  results. 

For  the  static  analysis  discussed  above  a uniform  positive  1.0  psi 
pressure  was  applied  to  the  SRM  model  throughout  the  combustion  cavity. 
When  the  first  longitudinal  acoustic  mode  is  active  in  the  combustion 
cavity,  the  pressure  mode  shape  calls  for  a positive  pressure  in  the  head 
end  at  the  same  time  the  pressure  is  negative  in  the  aft  end.  Thus, 
because  of  the  pressure  mode  shape,  a net  positive  thrust  can  be  produced 
by  the  aft  end  pressure  acting  over  the  aft  dome  area,  (less  the  nozzle 
opening),  and  acting  simultaneously  with  the  positive  thrust  caused  by 
the  integral  of  the  pressure  over  the  forward  dome. 

To  estimate  the  forward  dome  thrust  load,  the  1.0  psi  pressure  is 
multiplied  by  the  motor  cross  sectional  area: 


Ffwd  = (1.0)  7T  ( 72 ) 2 = 16286  lbs 

The  net  thrust  as  determined  above  for  a uniform  pressure  load  was: 


F fwd 


aft 


3570  lbs 


Therefore,  F , = 12716  lbs 

att 

Adding  Ff^j  to  Faft  to  account  for  the  pressure  mode  shape  gives  a total 
net  thrust  of 

Fn  = 16286  + 12716  = 29,002  lbs 

The  pressure  mode  shape  is  assumed  to  have  an  amplitude  of  + 1.0 
psi  in  the  forward  end  and  a -1.0  psi  in  the  aft  end. 
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The  total  thrust  of  29,000  lbs  acting  on  the  mass  of  the  SRB  would 
cause  a rigid  body  acceleration  of  (29,000  lbs/W),  where  W is  the  SRB 
weight.  If  the  SRB  weighs  1.26  x 10^  lbs,  then  a 15.25  Hz  oscillation 
would  occur  with  acceleration  amplitude  of 


29  x IQ-1 

3 = 1.26  x 106 


= 23  x 10‘3 


For  harmonic  motion  the  corresponding  displacement  would  be: 

y - -5-  - 23  x IQ"3  x 386,4  ,Q  . n-4  . 

X “ C°2)  (2  7T 15. 24)2  9,69  x 10  ln* 

The  above  results  show  that  the  SRB  could  be  expected  to  move  back  and 
forth  as  a rigid  body  at  15.25  cycles  per  second  with  a displacement  ampli- 
tude of  9.7  x 10'^  inches  in  response  to  the  first  longitudinal  acoustic 
mode . 

In  operation,  the  SRB  is  attached  to  the  ET  in  the  axial  direction  at 
a single  point  (node  303).  If  a rigid  SRB  were  attached  to  a rigid  ground, 
the  total  net  oscillatory  thrust  of  29,000  lbs  could  be  transmitted  through 
the  attach  point.  At  the  other  extreme,  if  the  ET  attach  point  were  very 
soft,  the  SRB  would  oscillate  near  the  rigid  body  displacement  of  9.7  x 10~~ 
inches  and  attach  point  forces  would  be  small.  Based  on  this  simplified 
rigid  body  analysis,  attach  point  forces  in  the  range  of  0 to  29,000  lbs 
would  be  expected.  The  corresponding  expected  displacements  would  be  in 
the  range  from  9.7  x 10“^  to  0.0  inches  with  the  larger  forces  correspond- 
ing with  the  smaller  displacements. 


To  obtain  a better  estimate  for  the  upper  limit  force  that  might  be 
expected,  consider  the  receptance  matrix  for  the  combined  ET  and  Orbiter, 
^RSS  • F°r  a crude  estimate  of  the  axial  force,  assume  that  only 
is  non  zero  in  the  equation: 


(Refer  to  equation  44  for  the  components  of  jUj  and  j F j . Then  the  first 
equation  from  the  above  matrix  equation  would  be: 


J303x  rHF303x 


The  matrices  [RrssJsyM  = RRSS  and  RRSS  ASYM  = RRSA  are  given  in 
Table  XIV.  From  the  table, 

(rll>SYM  = 7. 9596E-07  - 7.1221  E-08  i = 7.99  E-7  ^-5.1 


(ril)AsYM  = 4 . 6010E-07  - 4.1699  E-08  i = 4.62  E-7  ^-5.2 


Now,  suppose  the  maximum  free  body  displacement  of  9.7  x 10  inches  is 
applied  to  the  ET/Orbiter  combination  at  attach  node  303  in  the  axial  (X) 
direction: 


(F303X>SYM  = 77-)sym-  <U303X>  = jf 

1 9 

^OSX^ASWl  = (r  ,3  ^U303X^  = 2T 


9.7  x 10 
7.99  x 10 


-7  = 1214  lbs 


/t?  \ _ 1 /it  \ _ 9.7  x 10 

'•*  303X  ’aSTO  - (ru)ASYM  (U303X>  ~ 4.62  x 10'?  = 2x90  lbs 

For  comparison  with  the  numbers  calculated  in  these  rough  estimates, 
the  corresponding  values  from  the  computer  analysis  are  recalled: 


^U303X^SYM 

9.4  x 10 

(u303x  \s-m 

4.8  x 10' 

CF  ) = 

v 303X'SYM 

1306  lb 

^OSxVsYM 

1012  lb 

Based  on  the  relatively  good  agreement  between  the  rough  estimates  and  the 
computer  solutions,  we  can  justify  added  confidence  in  the  accuracy  and 
applicability  of  the  computer  solutions. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


This  analysis  program  has  resulted  in  estimates  for  the  attach  point 
forces  between  the  solid  rocket  motors  and  the  external  tank  and  has  there- 
fore met  the  program  objectives.  Good  agreement  between  hand  calculations 
and  computer  analysis  results  show  at  least  that  the  computer  results  are 
of  a reasonable  order  of  magnitude.  The  significance  of  the  maximum  attach 
point  force  of  1610  lbs  for  a + 1.0  psi  pressure  oscillation  level,  (or 
16,100  lbs  for  a + 10  psi  level)  must  be  determined  by  Space  Shuttle  engineers. 


Due  to  time  and  budget  limitations,  this  analysis  only  covered  the  first 
longitudinal  acoustic  mode  at  a zero  burn  time.  Future  work  could  include 
other  modes  and  other  burn  times.  Some  simplified  hand  calculations  could 
be  made  for  a transverse  mode.  The  desirability  of  additional  analysis 
work  could  best  be  determined  by  Shuttle  engineers  who  are  familiar  with 
space  shuttle  structure  capabilities  and  requirements. 
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